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TECHNICAL PAPER 


WELD STRESSES BEYOND ELASTIC LIMIT 
1. Material Discontinuity 

I. INTRODUCTION 


Many spacecraft structures are required to be verified with a minimum safety factor of 1.1 on 
yield and 1 .4 on ultimate strength [1] which, in the design phase, are both conservatively complied with 
in a linear fashion [2], As the structure progresses through the operational phase and prevailing environ- 
ments begin to exceed original design requirements, safety demands that current margins be affirmed. 
Determination of post-development safety factors on yield are linear and are usually accommodated, but 
ultimate strength predictions involve nonlinear considerations not commonly understood. The intent of 
this paper, the first in a series of three reports, is to provide mechanics insights to questions raised by Dr. 
George McDonough [3] on aluminum welded joint behavior beyond the elastic limit. 

The modified aft skirt structure of the solid rocket booster (SRB) was analyzed to predict safety 
margins using updated launch pad loads. Detailed global and substructural models were developed using 
a finite element method (FEM). The critical area was identified to be the aluminum weld joining the shell 
and a holddown post forging, because the weld material had a lower yield point and ultimate strength 
than adjacent parent metals. The weld was modeled with brick elements, five across the width, and the 
FEM code iterated the inelastic material properties in a piece-wise-linear routine [4], While the model 
predicted positive margins on the ultimate, the structural test article failed at the weld interface 8 percent 
below the required safety factor and at half the expected elongation. Twenty-four weld specimens were 
cut from the remaining skirt structure and uniaxially tested. All specimens were noted to fail at weld- 
parent interfaces. Average tensile design properties were reported and are included in Appendix A. 

Significant questions raised by this experience were: 

1 . What does consistent weld fracture at the interface imply? 

2. What is the effect of relative property differences between weld and parent materials on 
interface fracture? 

3. Is there a mechanics relationship between geometry and strength? 

4. What insights and considerations might be extended to FEM modeling? 

5. How much reliability was lost by the negative margin in post-development analysis? 

Of course there are many more questions, but these are vital for the revisit to the aft skirt design 
and for the choice of fabrication technology of future high performance structures. All of them have been 
addressed in some fashion in the past. Certainly, metallurgists have grappled extensively with the 
interface problem for years and are currently investigating it for the skirt failure. This paper responds to 
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the first four questions in some depth from a mechanics viewpoint and spills into a discussion on the last 
question. The most difficult task in relating the unknown interface phenomenon to a conceived model 
was not to find an elegant, or even an exact, approach, but to find any method at all that might describe 
the suggested physics. 


II. WELD INTERFACE STRESSES 


Perhaps the first suggestion that a discontinuity stress was acting at the weld interface came from 
Mr. Timothy P. Vaughn’s strain gage data [5] based on a uniaxial test of weld specimen (see Figure A-l 
in Appendix A). At 45-ksi stress, the gage on the weld center was indicating 0.022 strain, while gages at 
the plate and forging interfaces were reading about 4 percent higher strains. Parent plate strain away from 
the weld was less than 0.01. Metallurgical faults influencing consistent boundary fracture not 
excluded, it was thought that increased strain at both boundaries and fracture might be partially explained 
through mechanics. 

Referring to the stress-strain data (see Figure A-2 in Appendix A), at 45-ksi stress the plate is 
slightly strain hardened and the weld is far into it. The Poisson’s ratio for the plate is just a bit higher than 
elastic, but the weld is approaching 0.5. Consequently, under uniaxial tensile stress, lateral strains and 
related displacements at the interfaces should be sufficiently different to cause local discontinuity 
stresses. 

The approach was to model the test specimen, predict stresses and strains for the plate-weld 
boundary, and compare results with test gage data. Since discontinuity peak stresses dissipate over a very 
small span, a classical elasticity model was believed to define it more sharply than the FEM. Further- 
more, simplified assumptions could be more easily tracked and results intuitively qualified accordingly. 
The discontinuity model contrived is presented in Appendix B. It satisfies interface boundaries, and 
stress distributions are developed independently of material properties. Nonlinear properties required to 
determine displacements and strains are applied piece-wise linear. 


1. Uniaxial Specimen Stresses 

As a uniaxial tensile stress, cr x , is applied to the test specimen, the half thickness “h” contracts 
laterally a quantity v p along the plate and a greater quantity v w along the weld because of differences in 
Poisson’s ratio beyond the elastic limit (Fig. 1). 

This difference in lateral contraction may be expressed by 


Av = CT x hK^- 


a) 


where K is a simplified but effective relative material index defined by 
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It varies between 0 and -I, and resembles Figure 1 more as it approaches -1. Plate materials are 
referenced by a subscript (p) and weld properties are not subscripted. Moduli and Poisson’s ratios are 
calculated from Appendix A. 



Figure 1 . Lateral contraction under uniaxial stress. 

At the interface, £ = 0, the differential displacement of plate and weld must match according to 


Av(a x ) + v(q) = 0 


(3) 


where the first term is the relative displacement caused by the applied axial stress, a x , as determined from 
equation (1). The second term is the local stretching required from a dummy unit loading, q, acting adja- 
cent to the interface to conform to the displacement of the first term. This second term was derived using 
stress-strain relationships, and discontinuity factors D y and D x were developed in Appendix B to give 


v(q) _ fyh 

q q 


— [D y - |xD x ] 
E 


(4) 


Substituting equations (B-7), (1), and (4) into (3) gives the desired ratio of applied axial stress to applied 
dummy loading, 
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( 5 ) 


D = — = --%Fr + (l-ft)2 — sin(2Tot) . 
,rc i q llK a J 


Equation (5) is corrected for Gibb’s phenomenon at the interface x = d and is expressed in parametric 
form. Convenient parametric identities are 

T = — , H = -p , X = — , Y = — , and a = mil . (6) 

L L L L 

Using identities (6) into equation (B-7), stress discontinuity factors in the weld region are expressed in 
the following parametric form: 

D xy = 4 2 sin(aT) sin(aX) [H sinh(aY) - Y cosh(aY)] e' uM 

D y = T + 42 sin(aT) cos(aX) [(H + cosh(aY) - Y sinh(aY)] e uH 

D x = 4 2 sin(aT) cos(aX) [y sinh(aY) - (h - -2) cosh(aY)] e “" . (7) 


In applying these equations, note the two coordinate systems between Figures 1 and A-l . When £ = 0, 
x = d, and when £ — w/2, x = L. 

Finally, net distribution factors of stress components acting in the weld are related to the applied 
axial stress by dividing equation (7) by (5) and superimposing the applied stress to the axial equation: 


S 


xy 



D 


xy 


D 


<r.q 




1 + 



( 8 ) 


Superscript d denotes that discontinuity stresses are included and a x is the applied axial stress, as defined 
before. Using equation (8) and properties of weld test specimen at 45-ksi applied stress, 
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the computed distribution factors of stress components of equation (8) are graphically presented in 
Figure 2. 

Though these results cannot be adequately verified, there were two gage data points on the 
specimen surface that were thought to be correlatable; the weld midwidth strain gage and the interface 
gage which read 4 percent higher than the midwidth. Surface axial stress factor (S x ) predictions of 
equation (8) were obtained from the bottom chart of Figure 2 for y/h = 1 .0. The two test strain gage 
locations were fj/w = 0 and 0.5 corresponding to an interface estimated stress factor of 1 .07 and midweld 
of 0.96. 



Figure 2. Discontinuity factors of stress components in uniaxial specimen. 
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It was interesting to note that the midweld stress predicted at the surface was 4 percent less than 
the applied axial stress while subsurface stresses were equal to the applied stress. The net outcome was 
that the interface stress on the surface was a net 1 1 percent higher than at the weld midwidth. Referring to 
equation (A-3), an 1 1 percent variation of stress on the weld material translates to 

8e = — 5(11) = 48% strain (10) 

0.23 


increase at the interface. However, this strain was acting over a very narrow band of A£/w = 0.05, 
which was approximately 0.016 inches. The strain gage length was 0.125 inches which spanned over a 
weld strain of 0.024, and an interface strain of 0.036. Since changes in strain gage resistance in the two 
regions were in series, the interface strain gage net prediction was determined from the rule of mixtures 
to be 0.025, which was 4.5 percent higher than the midwidth location. 

This correlation was not entirely conclusive because it depended on dubious strain gage locations. 
If the gage span included a portion of the plate, the gage pickup would be less. If the plate material at the 
interface was softened by the welding heat, the gage reading would be higher, but the discontinuity factor 
would be less. The only point made clear so far was that, if an abrupt change of metallurgical properties 
occurs at the weld interface, then a stress riser could be induced and its characteristics would be 
consistent from weld to weld. This point was pursued further for other mitigating implications. 


2. Interface Discontinuity Stress Assessment 

If discontinuity stresses are inherent at weld interfaces, then welds should consistently fail at 
interfaces with predictable strength, and uniaxial tension test data should reflect the design property of 
this weakest link. Beyond that, designers must also know how discontinuity stress factors vary with other 
weld parameters. There is a need to know (1) how distortion energy theory of failure of combined 
stresses at the interface compares with the surface axial stress, and (2) how sensitive is the stress factor to 
material pairing using some relative material index. 

Distortion energy stresses are similarly derived in parametric form from the combined stress 
components of equation (8). Starting with the stress tensor. 


(S x -\) Sxy 
Sxy (S y — X.) 


= 0 , 


principal stress factors are determined, 


S,, = ± —*[(S X + Sy) 2 ~ 4(S X Sy — S x y)] 0 5 


(ID 
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These principal factors are substituted into the distortion energy theory 


S ( 2 ) = S] + S; - S,S 2 , 

so that 

S 0 = [S 2 + S 2 y - S,S 2 + 3S 2 y ] 05 . (12) 


The theory hypothesizes that a material will fail when the energy of distortion of combined stresses 
achieves the same energy for failure in simple tension. 

The distortion energy stress factor distribution over the weld is plotted in Figure 3. Again, as for 
the discontinuity factor in Figure 2, the surface energy factor away from the interface is less than the 
applied axial stress. This implies that test stress-strain relations based on surface strains are not conserva- 
tive. Those below the surface are reasonably close to the applied stress, as they should be, to comply with 
the distortion energy theory. The peak value occurs on the surface at the interface, which is only 1.5 
percent higher than the applied stress. 



Actually, distortion failure factors should be calculated at the failed applied stress of 49 ksi. 
Equation (1) shows the relative displacement between plate and weld to be directly proportional to the 
material property K. Since discontinuity stress factors are also directly related, they should increase with 
increase in K, as noted in Figure 4. The recalculated properties at weld failure of 49 ksi, index K = -0.5, 
and jlc = 0.5, produced an energy failure factor of S 0 = 1 .03 and a strain increase of 1 percent. Even for 
the worst case pairing of materials, the distortion failure stress factor at the interface of a uniaxial speci- 
men is not very sensitive to materials. 
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Figure 4. Material properties effects on uniaxial failure factor. 


3. Discontinuity Stresses in Plane Strain 

Consider a specimen having a very large width such as to ensure plane strain condition. How does 
the distortion failure factor vary in a different stress field? The only required modifications to the analysis 
are the relative material index, 

K ‘= + ’ (l3) 
and the ratio of applied axial-to-dummy stresses, 

D <r.q = if = + (1 -2\x) 2 sin(2Ta)J . (14) 


Incorporating these revised expressions to the above equations, the results are shown in Figure 5. 
Obviously, externally-imposed biaxial stress conditions have a decisive effect on the failure factor as it 
approaches ultimate stress and p = 0.5. This is a particularly meaningful discovery because the SRB aft 
skirt weld is in a biaxial stress field and the allowable weld strength was based on uniaxial test data. 
Consequently, the skirt failure is expected to always occur at the weld interface and at a lower strength 
than current predictions; perhaps 10 percent. 


8 




4. Implications to Other Fracture Issues 

There are four types of discontinuity stresses caused by abrupt changes in materials, geometry, 
loads, and temperature. Each may be solved piece-wise linear in the inelastic range, but whether any 
combination of the four occurring simultaneously may be linearly added needs to be established in fact. 
The aft skirt welds had three types of discontinuities at the holddown post fillet which compounded the 
failure predictability. 

If the same ultimate factor of safety is to be imposed on structures during development and post- 
operations, then more exploratory-type test programs in the inelastic range are required to provide 
appropriate design data and to verify nonlinear math models and prediction techniques. 

There are many conditions and opportunities for introducing metallurgical defects and distortions 
in aluminum weld interfaces [6] and particularly in large structures. Fortuitously, the weld peak failure 
stress occurs on the surface at the interface which can be inspected for flaws and other contributing 
deficiencies. However, weld discontinuity stresses do not seem to present a risk to fracture mechanics 
because the bandwidth of the peak stress field at the interface is much narrower than a practical flaw size. 
On the other hand, if the discontinuity failure factor is large, the fatigue life is reduced accordingly. 


III. SUMMARY 


When stress-strain relationships between weld and parent materials are appreciably different, 
discontinuity stresses at interfaces might be considered. If the weld consistently fails at the interface with 
no obvious metallurgical defect, then discontinuity stresses are a strong suspect. There are no straight- 
forward tests to detect them, but mathematical predictions followed by tests might provide implicit 
results and foundation to resolve this intriguing phenomenon. 
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The analysis presented and possible FEM techniques have their limitations. As stated before, the 
purpose for selecting the classical mechanics model was so discontinuity values and bandwidths might be 
more sharply defined. Perhaps iterating on nonlinear material properties might have improved resulting 
stress factors, but the motive for this analysis was to learn if there was a mechanics problem in the first 
place, and then to obtain insights on its implication in other related issues. 

Nevertheless, it may be concluded that 2219 aluminum weld material is sufficiently dissimilar 
from the parent metal such that discontinuity stress risers at the interface are inevitable, but its propensity 
to fracture is more acutely dependent on the externally-induced stress field. It is weak in uniaxial tension 
and significant in biaxial stress. Therefore, high performance welded shells (biaxial state) will always fail 
at weld interfaces. It is mildly sensitive to relative material properties of weld and parent materials. 

Though weld stresses vary throughout the uniaxial tensile test specimen, the ultimate applied 
stress should be indicative of the interface uniaxial strength. Also, strains vary along the weld specimen 
surface [5], slanting stress-strain data to the gage length and location. Therefore, it follows that weld 
strain data obtained from surface-mounted instrumentation on structural test articles may not be directly 
correlatable to FEM predictions beyond elastic range. 

The peak discontinuity stress occurs on the surface and at weld interfaces. If the peak stress riser 
in a biaxial field is large, the fatigue life is diminished accordingly. It has no effects on fracture 
mechanics because its bandwidth is smaller than most detectable flaws. 

This concludes responses to the first two questions raised in the introduction. 
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APPENDIX A 


WELD-PLATE MATERIAL PROPERTIES 


The 22I9-T87 aluminum structural properties of the plate and forging used in this document were 
obtained from the MIL-HDBK-C5. Weld properties were developed from 24 specimens excised from an 
SRB aft skirt which was constructed of 2219 aluminum plate welded to 2219 aluminum forging and heat 
treated to T87. The TIG weld had a thickness-to-weld ratio of 4.6 as shown in Figure A-l. These 
specimens were instrumented with electrical strain gages with 0.125-inch length. All specimens were 
uniaxially stressed in tension to failure in which 21 failed at the forging-weld interface while others failed 
at the plate interface. Average test results [5] of the weld and the C5 data are plotted in Figure A-2. 



Figure A-l. Weld specimen. 
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Figure A-2. Uniaxial test specimen properties. 
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Aluminum stress-strain relationship beyond the elastic limit may be represented by 


a = ae b 


(A-l) 


where “a” is the strength coefficient and “b” is the strain-hardening exponent. Both parameters were 
estimated from Figure A-2 plots and are noted as 


a b 

plate 125 0.24 

weld 106 0.23 (A-2) 


Using equation (A-l), the stress-strain variation relationship is calculated from 

42 = abe 1 -' ^2 , 

a e 


and letting Acr/cr = ba and A e/e = be gives, 


8a = b be 


(A-3) 


It may be reasonable to assume that since strain-hardening alters the stress-strain relationship according 
to equation (A-l), the same process should adjust the Poisson’s ratio between yield and plastic flow. A 
direct approach is suggested to let the Poisson’s ratio vary with the stress beyond the elastic range as 

B = By + tP'U By] f ~~~ °_ Y 1 • (A-4) 

[_au - a Y J 

where subscripts “U” and “Y” refer to ultimate and yield, respectively. Plate and weld Poisson’s ratio 
were assumed to be, 


|x Y = 0.3 and |Au = 0.5 , 


and related stresses in ksi are 
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a Y 

O’U 

plate 

32 

58 

weld 

26 

49 


Respective properties are applied to equation (A-4) to give the following desired relationship: 


|A p = 0.054 + 0.0077a 
|Xw — 0.074 + 0.0087a 


(A-5) 


The strain hardening effects on Poisson’s ratio may be related to strain by substituting the independent 
stress variable with equation (A-l). It has the form 

fx = C i + C 2 ae b 

The secant modulus is used as a relative stiffness parameter, 

E “f =[-^] l,b • (A ' 6) 

Using plate and weld properties given above, moduli are expressed respectively by 



Weld data presented is based on strain gages installed on weld surfaces at the center between plate and 
forging interfaces. 
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APPENDIX B 


DISCONTINUITY STRESS MODEL 


The stress problem is defined here as a welded plate which is symmetrical about the midwidth, 
w/2, and about the midthickness, h. The only external load acting on the weld is an axial tensile stress, 
cr x . The weld thickness shrinks at the center as the axial load increases, but it is partially constrained at 
the plate interface which imposes discontinuity stresses. 

The weld interface constraint produces a shear which is zero at the top and bottom surfaces, 
maximum just below the surfaces, and minimum at the midthickness. This condition may be approxi- 
mated by adding a dummy piece, as shown in Figure B-l , which allows the weld to be stretched through 
a dummy unit loading “q” over a length “d.” Since this dummy is not part of the weld structure, discon- 
tinuity stress distributions for x < d are not relevant to the analysis. 



x 


Figure B-l. Weld interface model. 

Figure B- 1 describes a simply connected contour with in-plane forces which makes it independent 
of material constants. The geometry and the discontinuous dummy loading “q” is symmetrical about the 
x and y axes. This is a classical problem [7], applicable to the generalized plane stress method. It consists 
of integrating the compatibility equation, 


2 d 4( E + 

dx 4 d\ 2 y 2 dy 4 


(B-l) 


where Airy’s stress function, <J>, is defined by three components derived from equilibrium: 


Ir'InjE^liyrk^ta aVU a t ; ii 
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(B-2) 




d 2 <\> 

'i 1 9 

dy~ 



•xy 


d 2 4> 

dxdy 


Equation (B-l) is of the form <l>(x,y) = XY whose solution is well known to be 


^*( x »y) — (C, cos Xx + C 2 sin Xx) (C 3 cosh Xy + C 4 sinh Xy + C$y cosh Xy + C^y sinh Xy) , 

(B-3) 


where C,...C 6 are arbitrary constants and X = rrnr/L for m = 1,2,... . Substituting equation (B-3) into 
(B-2), satisfying stress boundary conditions 


T xy = - - =0 at x = ± h and y = ± L , 

dxdy 

and allowing that sin Xx vanishes for a symmetrical loading of 

Q = A + I B cos 2= (B-4) 

1 L 

at y = ±h, then constants are determined to be 
C 1 = 1 and C2 = C4 = C5 = 0 » 

and stress components of equation (2) become: 


— X? cos Xx 

<r* ~ = 2Q Zj — — . t , |(Xh cosh Xh - sinh Xh) cosh Xy - (sinh Xh) Xy sinh Xy] 

dy- 2Xh + sinh 2Xh 


_ ^ 2( t> = 20 2 cos Xx 


°y — — 7 —2 Q Zj [(sinh Xh) Xy sinh Xy - (Xh cosh Xh + sinh Xh) cosh Xy] 

3x 2 2Xh + sinh 2Xh 3 2 2 


* xy 


a z <j> 

dxdy 


= 2Q 2 


sin Xx 


2Xh + sinh 2Xh 


[(Xh cosh Xh) sinh Xy - (sinh Xh) Xy cosh Xy] . (B-5) 
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Dummy uniform tension loading terms, A and B of equation (B-4), are found in the usual way: 


J* A dx = / q dx , 
or A = qd/L; 


-u 


B 


2 f d r L 2q 

i~ J Q cos -Vx dx + / (0) = — p- sin Xd , 

L -o XL 


(B-6) 


Evaluating A for integers 0 < m < 1 in each stress component, A becomes 0 for cr x and T xy , and A/2 for 
o-y. Assuming that Xh is large, then stress discontinuity factors are reduced to: 


_ -A» 

D y == = — + 4 X sin Xd cos Xx — [(Xh+ 1) cosh Xy - Xy sinh Xy] 

D x = — = 4 X sin Xd cos Xx [Xy sinh Xy - (Xh— 1 ) cosh Xy] 

XL 

D xv = — = 4 X sin Xd sin Xx ^ [Xh sinh Xy - Xy cosh Xy] . 

xy q xl 


External loading factor derived from equation (B-4) and discontinuity factors resulting from equation 
(B-7) are plotted in Figures B-2 and B-3 for m = 8. 
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Figure B-3. Discontinuity stress component factors. 
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